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Abstract. We study the production of cosmogenic neutrinos and photons during the extragalactic propagation of ultra-high- 
energy cosmic rays (UHECRs). For a wide range of models in cosmological evolution of source luminosity, composition and 
maximum energy we calculate the expected flux of cosmogenic secondaries by normalizing our cosmic ray output to experi- 
mental spectra and comparing the diffuse flux of GeV-TeV gamma-rays to the experimental one measured by the Fermi satellite. 
Most of these models yield significant neutrino fluxes for current experiments like IceCube or Pierre Auger. Furthermore, we 
discuss the possibilities of signing the presence of UHE proton sources either within or outside the cosmic ray horizon using 
neutrinos or photons observations even if the cosmic ray composition becomes heavier at the highest energies. We discuss the 
possible constraints that could be brought on the UHECR origin from the difl'erent messengers and energy ranges. 
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1. Introduction 

The origin of ultra-high-energy cosmic rays is one of the 
main topics in high-energy astrophysics. Their origin still re- 
mains unknown after decades of experimental efforts. Recent 
experiments such as the High Resolution fly's eye (HiRes, 
Abbasi et al., 2004), the Akeno Giant Air Shower Array 
(AGASA, Nagano et al., 1992) and primarily the Pierre Auger 
Observatory (Abraham et al., 2004) have increased the avail- 
able cosmic ray statistics above 10'^ eV, allowing the first solid 
studies from the point of view of the spectrum, the composi- 
tion and the arrival directions of ultra-high-energy cosmic rays 
(UHECRs). A sharp decrease of the UHECR flux above 3- 
4x10'^ eV possibly related to the theoretically expected GZK 
cut-off (Greisen, 1966; Zatsepin & Kuzmin, 1966) seems to be 
firmly established, whereas the interpretation of composition 
or arrival direction data as well as the consistency of the results 
reported by the different experiments are matters of intense de- 
bates. 

A multimessenger approach of the UHECR question is a 
promising way for solving the mystery of their origin. Indeed, 
shortly after the prediction of a cut-off at the highest ener- 
gies caused by the interactions of UHECR protons and nuclei 
with the cosmic microwave background (CMB) photons, the 
unavoidable associated production of UHE neutrinos through 
the decay of charged pions was pointed out by Berezinsky 
and Zatsepin (1969). Because these cosmogenic neutrinos can 
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travel from their production site without undergoing interac- 
tions or deflection, they were, soon after this pioneering work, 
considered as potentially interesting probes of the highest en- 
ergy phenomena in the universe. Their expected flux on earth 
was intensively calculated in the past decades for different 
astrophysical scenarios on the cosmological evolution of the 
sources, the composition or the maximum energy at the sources 
(see for instance Stecker ,1979; Hill & Schramm,1985; Engel 
et al., 2001; Kalashev et al., 2002; Seeker & Stanev, 2005; 
Hooper et al., 2005; Ave et al., 2005; Stanev et al., 2006; Aflard 
et al., 2006; Anchordoqui et al., 2007; Takami et al., 2009; 
Berezinsky, 2009; Ahlers et al., 2009; Kotera et al., 2010). 

Cosmogenic gamma-rays are also produced during the 
propagation of UHECRs. Unlike neutrinos, these very high- 
energy gamma-rays interact rapidly and produce electromag- 
netic cascades. As a result the universe is opaque to gamma- 
rays from a few hundreds of GeV to a few 10'** eV. Above 
10'^ eV the universe becomes more and more transparent to 
photons and very high-energy gamma-rays can propagate a 
few tens of megaparsecs without losing a great amount of en- 
ergy. As a result, these very high-energy cosmogenic gamma- 
rays were discussed in the literature either as signatures of 
the so-called Top-Down models (see for instance Protheroe & 
Johnson, 1996; Lee, 1998; Sigl et al., 1999; Semikoz & Sigl., 
2004) or as probes of UHECR acceleration in the local universe 
(e.g, Yoshida & Teshima, 1993; Protheroe and Johnson, 1996; 
Lee, 1998; Semikoz & Sigl., 2004; Gelmini et al., 2007a and 
2007b, Taylor & Aharonian, 2009; Taylor et al., 2009; Kuempel 
et al., 2009; Hooper et al., 201 1; Ahlers & Salvado, 201 1). 
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Because electromagnetic cascades, piling up below 100 
GeV, are produced during UHECR propagation and are likely 
to be associated with the production of cosmogenic neutri- 
nos, it has soon been realized that measurements of the diffuse 
gamma-ray background could allow one to put constraints on 
the cosmological evolution of the UHECR luminosity (Strong 
et al., 1973, 1974) and on the maximum allowable cosmo- 
genic neutrino fluxes (Berezinsky & Smirnov, 1975). Modern 
versions of this calculation were attempted using the EGRET 
measurements (Sreekumar et al., 1998; Strong et al., 2004) 
in Kalashev et al., (2002), Semikoz & Sigl (2004), Kalashev 
et al. (2007). More recently, the Fermi satellite measurements 
(Abdo et al., 2010) reported a gamma-ray background between 
100 MeV and 100 GeV lower than previously estimated us- 
ing EGRET data. The additional constraints allowed by this 
new measurement were discussed by Berezinsky et al. (2010), 
Ahlers et al. (2010) and more recently by Wang et al. (2011). 
These studies agree that the constraints implied by this new 
estimate of the gamma-ray background are much more strin- 
gent than those using EGRET data, but their conclusions differ 
on the impact of these new constraints on the observability of 
UHE neutrino fluxes (see below). In this paper we study the 
production of cosmogenic photons and neutrinos assuming dif- 
ferent astrophysical scenarios for the cosmological evolution 
of the source luminosity, the composition and different scenar- 
ios for the transition from galactic to extragalactic cosmic rays 
in the same spirit as in Kotera et al. (2010). We discuss the 
constraints brought by the Fermi measurements of the diffuse 
extragalactic background on the expectations for cosmogenic 
neutrinos or UHE photons detection. Furthermore, we discuss 
the possible neutrino or photon signatures from UHECR pro- 
ton accelerators either in the local universe or at cosmological 
distances (diffuse or point sources) even if the UHECR com- 
position becomes heavier, which is one of the plausible inter- 
pretations of the recent measurement of the longitudinal devel- 
opment of air showers (hereafter X„ax measurements) of the 
Pierre Auger Observatory (Abraham et al. 2010a). In the next 
section, we briefly introduce our code that simulates the prop- 
agation of UHECR protons and nuclei and our Monte Carlo 
tool for the development of intergalactic electromagnetic cas- 
cades. In section 3 we then calculate cosmogenic photons and 
neutrino fluxes by normalizing our simulations with the exper- 
imental spectrum measured by the Pierre Auger Observatory 
(Abraham et al., 2010b). We finally discuss these results and 
the possible constraints on the UHECR origin that can poten- 
tially be achieved by observing either point sources or diffuse 
contributions of cosmogenic secondaries. 

2. Modeling the propagation of cosmic rays and 
electromagnetic cascades 

In the following we consider the propagation of UHECR pro- 
ton and nuclei in the intergalactic medium, the production of 
secondary photons, pairs and neutrinos and the development of 
electromagnetic cascade within our Monte Carlo framework. 
In addition to the CMB photons, we consider in these calcula- 
tions all other main photon targets in the intergalactic medium. 
The universal radio background (URB) is considered in this 



study only for the development of electromagnetic cascades; 
for this purpose we conservatively use the spectrum given in 
Clark et al. (1970) and the cosmological evolution estimated 
in Lee (1998). Several authors have already treated the impact 
of different estimates of the radio background on the output 
of electromagnetic cascades (see for instance Gelmini et al., 
2008; Taylor and Aharonian, 2009). None of the conclusions 
of our studies should be significantly affected by the uncertain- 
ties related to the URB. In the infra-red, optical and ultra-violet 
range (hereafter IR/Opt/UV) we use the model of the density 
and cosmological evolution of the background photons pub- 
lished by Kneiske et al. (2004). We use the latest available up- 
date of their best-fit model, which proved to be compatible with 
the Fermi observations of distant blazars and GRBs (Abdo et 
al., 2010). For the sake of comparison, we also refer to the es- 
timate of Stecker et al. (2006) used in our previous works. This 
estimate was, however, challenged by the recent Fermi obser- 
vations suggesting a fainter photon density in optical and UV 
ranges. 

For the propagation of UHECR nuclei we used the code 
described in more detail in Allard et al. (2005,2006), where in- 
teractions of protons and nuclei were modeled using a Monte 
Carlo technique. All the relevant interaction and energy-loss 
channels were taken into account. Adiabatic losses and pair 
production were treated as continuous processes. Regarding 
the pair production we followed the cross section and inelas- 
ticities given in Rachen (1996), which we also employ for the 
scaling of these quantities with the mass and charge of the nu- 
cleus. The spectra of the produced pairs on the different photon 
backgrounds, which are in most cases the strongest contribu- 
tion to electromagnetic cascades (see below) were calculated 
following the work of Kelner and Aharonian (2005). The nu- 
cleon pion production, which leads to the production of pho- 
tons, neutrinos and electrons or positrons, was simulated using 
the SOPHIA event generator (Mucke et al., 2001). In the case 
of nuclei, the different channels of the giant dipole resonance 
(GDR) were modeled using the theoretical calculations from 
Khan et al. (2005) and parametrizations from Rachen (1996) 
for nuclei with mass A< 9. Cross sections and nucleon yield for 
the quasi-deuteron and pion production (baryonic resonances) 
were also treated using Rachen (1996). We also followed his 
treatment for the nuclear transparency and the kinematics of 
the produced pions. Our Monte Carlo code allows us to propa- 
gate protons or nuclei from their source to the observer, follow- 
ing the evolution of the energy, mass and charge and the pro- 
duction of secondary photons, pairs and neutrinos. Neutrinos 
are supposed to suffer only from adiabatic losses during their 
propagation, while e^e" pairs and photons eventually initiate 
electromagnetic cascades. 

The latter were calculated using a Monte Carlo code com- 
puting the most relevant interaction channels (i) for photons: 
the e^e" pair production (Peskin & Schroeder, 1995) and the 
double pair production (following Lee, 1998), (ii) for electrons 
and positrons: the inverse Compton scattering (Jones, 1968; 
Zdziarski & Svensson, 1989) and the triplet pair production 
(Haug, 1975; Mastichiadis, 1991; Anguelov et al., 1999). The 
competition between these processes was treated stochastically, 
whereas the synchrotron and adiabatic losses were treated as 
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continuous processes. Cascades were calculated in a thin grid 
of electrons and photons with a primary energy between 100 
MeV and 10^' eV (with a 0.05 logarithmic step) and a grid 
of redshift between 10"^ and 8 (with a 0.05 logarithmic step). 
Five-hundred cascades were calculated for each combination. 
To limit excessive computation times, we applied a simple 
thinning algorithm as soon as the cascade particles fell below 
0.01% of the primary energy. When a secondary photon, elec- 
tron or positron is produced in our UHECR propagation Monte 
Carlo code, a pre-calculated electromagnetic cascade is ran- 
domly chosen at the corresponding energy and redshift and the 
electromagnetic particles reaching the earth are stored along 
with cosmic rays and secondary neutrinos. We defined the ef- 
fective loss length as the mean distance needed for the most en- 
ergetic particle of the cascade (which is most of the time not the 
same as the primary particle that initiated the cascade) to fall 
below 37% of the energy of the primary particle. The effective 
loss length was computed for photon-initiated cascades using 
our Monte Carlo, the result is displayed in Fig. [T] Although 
numerically different from the loss length calculated in Lee 
(1998) because of the different definition used, the same fea- 
tures in the energy evolution are reproduced. The successive 
roles played by the different photon backgrounds in the opacity 
of the universe to electromagnetic particles are visible. In the 
following we assume an extragalactic magnetic field of 10"^ 
nG in all our calculations. Although we do not discuss the in- 
fluence of the field strength in more detail, it was previously 
shown that fields below 1 nG do not significantly modify the 
expected photon flux above 10'^ eV (see for instance Gelmini 
et al., 2008) and we checked that the expected GeV-TeV diffuse 
flux was not strongly modified by magnetic fields of these mag- 
nitudes. We keep in mind, however, that fields stronger than ~1 
nG might strongly modify UHE photon fluxes as soon as UHE 
cascades start developing. 

3. Principle of the calculation 

In the next sections we calculate secondary photons and neu- 
trino fluxes for different astrophysical models in terms of com- 
position, maximum acceleration energy at the sources, and cos- 
mological evolution of the cosmic ray luminosity. The cos- 
mic ray spectrum measured by the Pierre Auger Observatory 
(Abraham et al., 2010) allows us to normalize our UHECR and 
secondary neutrino and photon fluxes and to choose spectral 
indices compatible with the overall shape of the UHECR spec- 
trum. In most cases, we will consider maximum energies at 
the source Emax = Z x 10^'^'^ eV (i.e., we assume that pro- 
tons and nuclei are produced above the pion production thresh- 
old) and perform an exponential cut-off of the injection above 
this energy. Sources are continuously distributed between a 
redshift z = 8 and a minimum distance of Dn,in - 4 Mpc. 
We deliberately did not perform a complete statistical analy- 
sis to choose these spectral indices, as was done in Ahlers et 
al. (2010). Indeed, the compatibility of the different models 
with respect to experimental spectra can slightly differ when 
one compares them to the Auger or HiRes data. For instance in 
the framework of the pair production dip model (Berezinsky et 
al., 2006) good fits to experimental data can be found more eas- 
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Fig. 1. Effective loss length (i.e., distance needed for the most 
energetic particle of the electromagnetic cascade to have less 
than 37% of the energy of the primary photon) as a function of 
the energy of the photon that initiated the cascade. The spread 
(1 RMS) is also displayed. 



ily with the HiRes estimate of the UHECR spectrum than with 
the Auger estimate. In this context, we simply use the spectral 
indices that yield the best compatibility with either the Auger 
or HiRes spectrum. A more detailed study of the compatibility 
of the different models to experimental data, where a detailed 
description of the systematic uncertainties of both experiments 
would be necessary, would not qualitatively change any of the 
conclusions in our study however. 

For our hypotheses on the cosmological evolution of the 
source luminosity we follow Kotera et al. (2010) and consider, 
besides a stationary case (hereafter labeled as "uniform"), the 
cosmological evolution of the source luminosity correspond- 
ing to the star-formation rate (SFR) estimated in Hopkins and 
Beacom (2006) as well as the evolution of powerful radio 
galaxies (labeled below FR-II) derived from Wall et al. (2005). 
Other hypotheses on the luminosity evolution were used in 
Kotera et al. (2010) such as alternative models of the star- 
formation rate or different estimates of the luminosity evolution 
of gamma-ray bursts; these cases were yet found to give sec- 
ondary fluxes very similar to those produced in the SFR case 
and therefore were not considered in the present study. 

We compared our numerical model of electromagnetic cas- 
cades with previous works by Ahlers et al. (2010) and Kalashev 
et al. (2009). The comparison is displayed in fig. l2] where the 
photon fluxes are normalized by choosing a similar normaliza- 
tion for the different curves with respect to the corresponding 
UHECR flux at lO'^ eV. To provide a consistent comparison 
with the two studies cited above, we used both Stecker et al. 
(2006) as well as Kneiske et al. (2004) estimates of IR/Opt/UV 
background photons. None of these backgrounds were used in 
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Fig. 2. Expected fluxes of secondary photons for a pure pro- 
ton composition, a spectral index jS - 2.0, a maximum en- 
ergy Emax = 10^' eV, compared to estimates by Kalashev et al. 
(2009) and Ahlers et al. (2010). The units are arbitrary but the 
different fluxes were consistently normalized with respect to 
their corresponding cosmic ray flux at 10'^ eV (these UHECR 
fluxes are omitted in this figure, but agree very well). The 
IR/Opt/UV backgrounds estimated by Kneiske et al. (2004) and 
Stecker et al. (2006) are used for these comparisons (see text). 

Ahlers et al. (2010), although their results appear to be very 
close to those we obtained using Kneiske et al. (2004). The 
estimate of Stecker et al. (2006) was used in Kalashev et al. 
(2009). In both cases the comparison agrees very well with pre- 
vious calculations. One can see the influence of the modeling 
of the IR/Opt/UV background on the expected photon flux in 
the Ge V-TeV region. The fainter optical and UV background of 
Kneiske et al. (2004) triggers a pile-up of secondary photons at 
higher energy than the background of Stecker et al. (2006). The 
modeling of the optical and UV backgrounds appears to have 
some influence on the shape of the resulting diffuse gamma-ray 
fluxes in the GeV-TeV region. Although it is important to stress 
that these photon backgrounds are not perfectly constrained 
at the moment, the estimate of Kneiske et al. (2004) shows a 
much better compatibility with the recent Fermi observations 
(Abdo et al., 2010). At the highest energies, the modeling of 
the IR/Opt/UV backgrounds becomes totally irrelevant for sec- 
ondary fluxes, our simulations agree well with previous works, 
especially with Ahlers et al. (2010). 

Before presenting the results of our study, we will compare 
the production of secondaries for different species of UHECR 
primaries, namely protons, helium, and iron nuclei (we will 
mainly concentrate on the production of photons because a 
detailed discussion about neutrinos can be found in Allard et 
al., 2006, and Kotera et al., 2010). A comparison between dif- 
ferent primary composition is displayed in Fig. l3] One can 
see the cosmic ray, neutrino and photon outputs for a source 
luminosity distribution following the SFR for three composi- 
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Fig. 3. Top; cosmic ray, neutrino (summed over all flavors) and 
photon spectra assuming three compositions: pure proton, pure 
helium and pure iron at the source, a source spectral index 
I3-2.Q and a maximum energy at the source Emax(Z) - Zx 10^"'^ 
eV The same cosmic ray luminosity between 10^^ eV and 
Emax(Z) is assumed. The contribution of the main fragment is 
shown in thin dashed lines. Bottom: same as top panel but as- 
suming a source spectral indeXjS=2.3. 



tions: pure protons, pure helium, and pure iron. The neutrino 
fluxes shown in all figures are summed over all flavors. We 
assume that the same luminosity is injected between lO'^ eV 
and E,nax(Z) = Z x lO^"'^ eV for the three composition mod- 
els. In the top panel a spectral index (3 - 2.0 is assumed. The 
resulting photon spectra appear very similar whatever the pri- 
mary composition assumed. Indeed, in the case of helium or 
iron primaries, most of the photon flux is caused by secondary 
nucleons. Above ~2-3xlO'^eV, secondary nucleons are 
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emitted within a few Mpc from the sources owing to photodis- 
integration of primary nuclei by CMB photons. Below the pho- 
todisintegration threshold with CMB photons, nuclei interact 
with far-infrared photons, although the photodisintegration is 
much slower (see the evolution of the mean free path in Allard 
et al., 2006), it remains very efficient on cosmological scales 
and most of the nuclei emitted by very distant sources release 
most of their nucleons during propagation. Assuming a com- 
plete photodisintegration, secondary nucleons emitted by par- 
ent nuclei in the energy range [AxE, Ax(E+dE)] are produced 
in the energy range [E, E-i-dE] . Their relative abundance with 
respect to primary nuclei in the same energy range is A^"^ and 
so is the relative total energy they carry. With a spectral index 
jS = 2.0 secondary nucleons are then as numerous as primary 
injected nuclei at the same energy, provided that the parent nu- 
cleus is totally photodisintegrated (which is almost always true 
for nuclei interacting with CMB photons unless the source dis- 
tance from the observer is below ~ 5 Mpc). Thus, they carry 
the same total amount of energy as primary protons would do, 
as soon as the same luminosity is assumed. At the highest ener- 
gies, the gap between the proton and nuclei cases is mostly due 
to the fact that we assumed a maximum energy proportional 
to the charge of the nucleus (to be consistent with the assump- 
tion that we make below) and would mostly disappear if we 
had assumed a scaling with A (note that in this case, the iron 
case would remain slightly lower owing to the higher absorp- 
tion probability for the produced pion). The contribution of the 
residual nucleus to the low-energy photon flux (shown in thin 
dashed lines) is quite low. This contribution is mainly caused 
by the pair production mechanism. For nuclei the loss length 
decreases as Z^/A at a given Lorentz factor. This process is 
thus quite efficient for iron nuclei above ~ 6 x lO'** eV (at z=0, 
see Allard et al., 2006). However, unlike in the proton case, this 
process has to compete with photodisintegration and is rapidly 
overwhelmed by GDR interactions with CMB photons (above 
~ 3 X 10^*^ eV for iron nuclei at z=0). As a result, it is only 
efficient on a short energy range. 

As soon as the spectral index becomes softer than 2.0 (see 
bottom panel of Fig.[3]where a spectral index/? = 2.3 is consid- 
ered), the injected luminosity decreases with energy, the pure 
proton and compound nuclei cases are no longer equivalent and 
the energy transferred to neutrinos and electromagnetic parti- 
cles becomes lower as the mass of the parent nucleus increases 
(according to the scaling in A^"^) . Because the threshold for 
the pair production of nuclei is lower than for GDR interactions 
with CMB photons, the relative contribution of this process 
(and thus the contribution of the parent residual nucleus) be- 
comes higher for softer spectral indices. Unless assuming very 
soft spectral indices, expected fluxes of secondaries remain in 
the same order as in the proton case. A similar discussion can 
be found in Ahlers and Salvado (20 11). 

As can be seen in Fig. [3] for the same assumed luminosity, 
the flux of cosmic ray expected on earth is higher in the case 
of a pure iron composition. This is because protons and then 
helium have a lower threshold for pair production interactions 
and as a result suffer faster energy losses than heavier nuclei 
(at least between ~ 10"* and 3 10^" eV see Aflard et al., 2006 
and 2008). Their flux is then more attenuated in this energy 



range. To reach a certain cosmic ray flux normalization for in- 
stance at 10^^ eV, the required luminosity is usually lower if 
a heavy composition is assumed, and then contributes to the 
usually lower secondary fluxes found when assuming a heavy 
composition (see below). 

4. Diffuse neutrino and photon fluxes for different 
astrophysical models 

In the same spirit as Kotera et al. (2010), we studied different 
models assuming various compositions and involving different 
interpretations of the transition from galactic to extragalactic 
cosmic rays. In all models in this section a maximum energy 
Emax = Z X 10^"'^ eV is assumed, i.e. we consider cosmic rays 
accelerated above the threshold of pion production with CMB 
photons. We injected cosmic rays from 10"' eV to Emax(Z) with 
various spectral indices (choosing for each source evolution hy- 
pothesis the spectral index that agrees best with experimental 
data) and assumed an exponential cut-oflF of the injection above 

^maxv-^/- 

4.1. dip model 

We started our study with the pure proton pair production dip 
model (Berezinsky et al., 2006). This model is very specific 
because here the ankle is not related to the transition from 
galactic to extragalactic cosmic rays, but is instead the sig- 
nature of proton interactions with CMB photons via the pair 
production mechanism. The resulting cosmic ray, neutrino and 
photon spectrum that we obtained for the three hypotheses on 
the cosmological evolution of the source luminosity are pre- 
sented in Fig. [4] Calculations are shown for both estimates of 
the IR/Opt/UV from Stecker et al. (2006) and Kneiske et al. 
(2004). Evidently, the photon (at low-energy) as well as the 
neutrino fluxes are very dependent on the assumed evolution 
of the source luminosity (see for instance Seckel and Stanev, 
2005). Indeed, both low-energy photons and neutrinos come 
from cosmological distances and can benefit from a lower in- 
teraction threshold and denser photon backgrounds at high red- 
shifts. In the case of a strong evolution of luminosity with red- 
shift, the sources with high redshift have a relatively greater 
weight, which immediately turns into higher secondary fluxes. 
This is not the case for high-energy photons, because these par- 
ticles always cascade to low energies unless they are produced 
in the local universe. The luminosity evolution is thus barely 
relevant for high-energy gamma-rays, but tends to slightly at- 
tenuate the production of high-energy photons because a lower 
weight is given to local sources. On the other hand, the spec- 
tral indices required to fit the experimental data are harder for 
strongly evolving luminosities with redshift, which means that 
relatively more luminosity is injected at high energies. In the 
dip model the two effects compensate almost exactly. Therefore 
the evolution of the source luminosity has basically no visible 
effect on the high-energy photon flux. 

When comparing the influence of the IR/Opt/UV model 
used, one can see that it mainly affects neutrino fluxes in the 
PeV region. This point was already discussed in Kotera et al. 
(2010): the fainter photon background in the optical and UV 
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Fig. 4. Cosmic ray (markers), neutrino (dashed lines) and pho- 
ton (solid lines) spectra (E^ x dN/dE) for the dip model com- 
pared to Auger spectrum (Abraham et al., 2010; open circles) 
and the Fermi diffuse gamma-ray spectrum (Abdo et al., 2010; 
black squares). The contribution of the pion mechanism to the 
photon spectrum is shown (dashed lines). The chosen spectral 
indices are/? = 2.6 for the uniform case (no evolution), 2.5 for 
SFR and 2.3 for FR-II. The results were computed assuming 
the IR/Opt/UV background estimate from Stecker et al., 2006 
(Top) and Kneiske et al., 2004 (Bottom). In the top panel the 
Auger 90% C.L integrated upper limit (2 years) for tau neutri- 
nos assuming a pure E"^ neutrino spectrum is also shown for 
comparison (Abraham et al., 201 1; the line represents the cen- 
tral value and was multiplied by 3 assuming a complete mixing 
of the neutrino flavors). The equivalent IceCube limit (IC-40, 
red thick-dashed line) is also shown (Abbasi et al., 20 11). 



range from the estimate of Kneiske et al. (2004) leads to neu- 
trino fluxes a factor of ~ 2 lower at 10'^ eV and dropping much 
faster below this energy. For both of the background models the 
expected low-energy photon fluxes significantly overshoot the 
diff'use photon flux measured by Fermi in the scenario of a FR- 
II evolution of sources. Constraints seem to be more stringent 
using the photon background by Kneiske et al. (2004), favored 
by the Fermi observations (Abdo et al., 2010) and in this case 
the photon flux in the SFR evolution case appears to be very 
close to the Fermi bounds. Here, we confirm previous results by 
Berezinsky et al. (2010) and Ahlers et al. (2010), claiming that 
in the framework of the dip model, the Fermi measurements of 
the diffuse gamma-ray flux actually involve strong limitations 
on the expected cosmogenic neutrino fluxes. By themselves, in- 
deed, ruling out basically all models that yields neutrino fluxes 
higher than the SFR model, they imply neutrino fluxes almost 
an order of magnitude lower than the upper limit of the Pierre 
Auger Observatory (see Abraham et al., 2009; Tiffenberg et 
al., 2009; Abreu et al., 2001 and Fig.Hl) and even lower than 
the current limits from the IceCube collaboration (Abbasi et 
al., 2011). Constraints obtained from the Fermi measurements 
can be somewhat dulled by invoking a low-energy cut mecha- 
nisnjjthat would leave the UHE neutrino flux unchanged while 
decreasing the pair production contribution (see below) to the 
diff'use gamma-ray flux. However, this would be at the expense 
of the cosmogenic neutrino flux between 1-100 PeV (see AUard 
et al., 2006). 

4.2. Mixed composition transition models 

We now consider the mixed composition model from Allard et 
al. (2005). There, the composition at the extragalactic sources 
is assumed to be similar to that of low-energy galactic cosmic 
rays. In this case, a pair production dip is no longer possible 
because of the significant contribution of nuclei to the source 
composition, and one can fit the cosmic ray spectrum down 
to the ankle (which is in this case the signature of the end of 
the transition from galactic to extragalactic cosmic rays) with 
harder spectral indices than for the dip model. Results are dis- 
played in the top panel of Fig.lS] One can see that in this case, 
as previously shown in Allard et al. (2006) and Kotera et al. 
(2010), the high-energy neutrino and UHE photon fluxes are 
very similar to the one obtained for the dip model. At PeV en- 
ergies, the neutrino fluxes are, however, much lower because of 
the harder spectral index required to fit the experimental data 
which leads to lower injected luminosities at low-energy. 

The constraints implied by the Fermi diff'use flux appear 
to be much less stringent for the mixed composition model 
than for the dip model. Only the FR-II source evolution model 
seems to be constrained by slightly overshooting the Fermi 
flux, while the low-energy photons produced in the SFR case 
are safely below the bounds. For the mixed composition model, 
the bounds given by Fermi are only constraining the most opti- 



' A change of the spectral index below ~ 10"* eV to a harder value 
owing to a change of the acceleration regime at the source that allows 
one to limit the luminosity injected at low-energy, see Berezinsky et 
al. (2006) 
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mistic neutrino flux expectation at the same level as the current 
upper limits from the Pierre Auger Observatory. The recent up- 
per limit of IceCube clearly rules out the specific FR-II sce- 
nario and already constrains less optimistic scenarios between 
the FR-II and the SFR case. We will return to this discussion 
below, but we can already point out that the cosmogenic neu- 
trino fluxes implied by the diffuse gamma-ray flux are strongly 
model-dependent and the possibility of observing high cosmo- 



genic neutrino fluxes is not excluded, confirming the conclu- 
sion of Ahlers et al. (2010). 

The bottom panel of Fig. |5] shows the case of a pure iron 
composition at the source. Although very different from the 
mixed composition case from the point of view of the source 
composition, the pure iron case has basically the same im- 
plications for the interpretation of the ankle. The pure iron 
case shows stronger variations of the UHE neutrino and pho- 
ton fluxes with the cosmological evolution assumed than other 
models. This is caused by a stronger variation of the spectral 
index required to fit experimental spectra. The uniform case 
shows a flux of UHE neutrinos and photons much lower than 
the pure proton case owing to the soft spectral index (2.4) re- 
quired (although it is harder than 2.6 for the dip model) and, as 
mentioned in the previous section, the lower luminosity needed 
for an iron composition (see below). For stronger evolution, 
harder spectral indices are needed and the gap between the pure 
proton and pure iron cases becomes smaller (a factor ~ 2 in the 
SFR case) and eventually the fluxes are similar for the FR-II 
evolution where, however, a /? = 1.6 spectral index is required. 
From the point of view of the Fermi, Pierre Auger Observatory 
and IceCube diffuse flux implications, the pure iron case ap- 
pears very similar to the mixed composition case. 



4.3. Classic ankle and late-transition models 



We also considered later transition models such as the classic 
ankle model considered for instance by Wibig and Wolfendale 
(2005). For that purpose a/3 = 2.0 spectral index is assumed for 
all evolution models. Results are displayed in the top panel of 
Fig.l6] The experimental spectrum can only be matched above 
10'^ eV for this model, and high-energy neutrinos expectations 
are then not higher than for the dip model (see Takami et al., 
2009; Kotera et al., 2010) despite the harder spectral, owing to 
the lower luminosity needed to match only the highest energy 
point of the spectrum (see discussion below). The constraints 
brought by the Fermi diffuse flux are then much looser than for 
the dip model and comparable to those of the mixed composi- 
tion model. A more stringent constraint comes from IceCube 
which is now able to rule out the highest neutrino fluxes ex- 
pected in the FR-II scenario. 

Finally, we studied a very late transition model, assuming 
an extremely hard spectral index y6 = 1 .0 for all evolution sce- 
narios (see the bottom panel of Fig. [6]). This injection spectrum 
indices are for instance expected for an acceleration by strong 
electric fields in the environment of young magnetars (Epstein 
et al., 1999; Arons, 2003). In this case, this type of extragalac- 
tic sources only significantly contributes to the highest energy 
end of the spectrum and the required luminosity is then very 
low, and UHE neutrino and photon outputs do not exceed the 
previous cases. But the Fermi diffuse flux does not constrain 
any of the evolution models by itself, whereas the neutrino 
fluxes obtained in the FR-II scenario are already constrained 
by IceCube. 
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Table 1. Cosmic -ray luminosity per unit volume (at z=0, 
in unit lO'^'^ergyr-' Mpc-^) above 10'^ eV (resp. above lO" 
eV) required for the different models displayed in Figs. 4]|6 
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Fig. 6. Same as Fig.|4| Top: a spectral index p - 2.0 is assumed 
for all source evolution scenarios to reproduce the classic an- 
kle transition model (see for instance Wibig and Wolfendale, 
2005). Bottom: a spectral index/? = 1.0 is assumed to produce 
a very late transition best motivated for magnetar-type sources 
(Epstein et al., 1999; Arons, 2003). 



5. Observability of diffuse fluxes of UHE neutrinos 
and photons 

5.1. Fermi constraints on cosmogenic neutrino fluxes 

In the previous section, cosmogenic secondary fluxes for sev- 
eral astrophysical models were calculated. We have shown that 
as previously pointed out by Ahlers et al. (2010), the impli- 
cations of the Fermi measurements of the diffuse gamma-ray 
flux are model-dependent. This model dependance can eas- 



Spectral indices of sources can be found in the corresponding 
figures. 



model 


uniform 


SFR 


FR-II 


pure proton dip 


11.8(47.4) 


7.5 (24.0) 


4.24(9.1) 


mixed comp. 


7.1(15.1) 


4.1 (6.1) 


2.4 (2.8) 


pure iron 


4.0 (9.4) 


1.8(2.2) 


1.6(1.8) 


classic ankle 


2.2 (2.9) 


2.2 (2.9) 


1.6(2.2) 


late transition 


1.1(1.1) 


0.7 (0.7) 


0.7 (0.7) 


mixed low E^ax 


5 (7.5) 


2.6 (2.6) 


2.2 (2.2) 



Table 2. Cascade energy densities Wcas (in unit eV cm ^) for 
the models displayed in Figs.|4]and|5] 



model 


uniform 


SFR 


FR-II 


pure proton dip 

mixed comp. 

pure iron 


1.3 X 10-' 
4.8 X 10-** 
9.5 X 10-' 


5.0 X 10-' 

1.7 X 10-' 

6.8 X 10-** 


1.7 X 10-« 
6.9 X 10-' 
5.1 X 10-' 



ily be understood by noticing that large parts of the GeV-TeV 
gamma-ray flux produced by UHECRs are generated by the 
pair production mechanism. The energy threshold of this pro- 
cess (for the interaction with CMB photons, which is respon- 
sible for most of the energy transfer between cosmic rays and 
secondaries) is much lower than the pion production mecha- 
nism that can either produce neutrinos, photons, or electron and 
positrons. In other words, the neutrino flux is related only to the 
subdominant pion-induced component of the diffuse gamma- 
ray flux, which means that similar UHE neutrino fluxes can 
be associated with a wide range of GeV-TeV gamma-ray flux 
depending on how high the flux produced by the pair produc- 
tion mechanism is. For a given normalization of the cosmic 
ray spectrum at 10'^ eV, the contribution of the pair production 
mechanism to the energy transfer to electromagnetic particles 
will ultimately depend on the cosmic ray luminosity injected 
at low-energy, which is related to the spectral index needed for 
a given model to reproduce experimental data. The very spe- 
cific interpretation of the ankle, as a consequence of the in- 
teraction of protons with CMB photons, requires soft spectral 
indices for the dip model with respect to other models con- 
sidered here and thus involves a higher luminosity injected at 
lower energy (i.e. below 10'^ eV). The cosmic ray luminosi- 
ties above 10'^ and 10"^ eV needed for the different models 
we considered are summarized in Table. [T] Evidently the dip, 
indeed, does require higher luminosities than other models, es- 
pecially in the energy range where the pair production mecha- 
nism is dominant (e.g above a few 10^^ eV depending on the 
assumed evolution of the source luminosity). These higher lu- 
minosities explain the larger diffuse gamma-ray flux at low- 
energy and the fact that the dip model is more constrained 
by the Fermi diffuse gamma-ray flux than the other models 
with a maximum allowed flux almost an order of magnitude 
lower, unless a low-energy-cut mechanism is assumed. These 
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previous results can be summarized concisely by comparing 
the cascade energy density Wcas (Berezinsky & Smirnov, 1975) 
for some of the models we discussed in the previous section 
(shown in Tab. l2]) with the maximum cascade energy density 
w™" = 5.8 X 10"^ eVcm-3 (Berezinsky et al., 2010) implied 
by the Fermi measurement of the diffuse gamma-ray flux. For 
models other than the dip model, which yield equivalent UHE 
neutrino flux expectations for a given evolution scenario, Fermi 
bounds are only saturated in the FR-II evolution scenario. This 
means that, there, the constraints obtained from neutrino ob- 
servations are already more stringent than those deduced from 
the observations of the diffuse flux of gamma-rays. Hence, the 
Pierre Auger Observatory and neutrino observatories will cer- 
tainly keep providing us with interesting constraints in the next 
few years. In particular, IceCube upper limits will probably 
reach the level of the predictions in the SFR source evolution 
scenario for most of the models we presented in the previous 
section within ~ 3-4 years. The Fermi constraints can be made 
more severe however, by identifying significant contributions 
of astrophysical sources to the diffuse gamma-ray flux. On this 
subject, the Fermi collaboration (Abdo et al., 2010) recently re- 
ported a contribution on the order of ~ 20% dominated by BL 
Lac objects. Current levels of this contribution do not strongly 
modify our conclusions yet. 

The discussion on the implications of the Fermi diffuse 
gamma-ray flux is different for cosmogenic neutrinos in the 1- 
100 PeV range. Indeed, in this energy range cosmogenic neutri- 
nos are mostly produced by interaction of cosmic ray between 
10'^ and 10'^ eV with higher energy background (i.e., infrared, 
optical or ultra-violet). Strong neutrino fluxes also require high 
luminosities in this energy range [i.e. soft spectral indices), then 
the predicted fluxes for the dip model are usually much higher 
than for the other models (see Kotera et al., 2010). The photon 
background density evolution with redshift in the IR/Opt/UV 
range is weaker than for CMB photons, and moreover, spectral 
indices required to fit the cosmic ray spectrum are becoming 
harder for scenarios of a strong evolution of the source lumi- 
nosity. As a result, the difference between the different source 
evolution scenarios in terms of neutrino flux in the 1-100 PeV 
range is smaller than for UHE neutrinos. Unlike for UHE neu- 
trinos, the difference between transition models is larger than 
the difference between source evolution scenario (e.g. the neu- 
trino flux for the dip model with SFR evolution is higher than 
for the mixed composition model and FR-II evolution). Then 
the highest allowed cosmogenic neutrino fluxes are those pre- 
dicted for the dip model and the SFR evolution. In this case, the 
event rate for IceCube was estimated in Kotera et al. (2010) be- 
tween 0.8 and 0.2 events per year using the Stecker et al. (2006) 
estimate of the IR/Opt/UV background. One would have to de- 
crease this rate by a factor ~2-3 if using Kneiske et al. (2004) 
results. 



5.2. Observability of UHE cosmogenic photons 

Because the cosmological evolution of the source luminosity 
does not strongly influence UHE cosmogenic photon fluxes 
(see above), the Fermi diffuse flux does not put any constraint 
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Fig. 7. Top: cosmic ray (solid lines), and photon (dashed lines) 
spectra (E- xdN/dE) compared to Auger spectrum for different 
models assuming a SFR source evolution : dip model, jS - 2.5, 
Eniax = 10-"-^ eV; late-transition model, /3 = 1.0, E,™x = lO^"'^ 
eV and 10^' eV (short dashed lines); mixed composition, fi - 
2.1, E„ax = 10^1 eV; pure iron, /3 = 2.0, En,ax = lO^"-^ eV 
Bottom: Integrated flux for the same models. Upper limits on 
the UHE photon flux from the Pierre Auger Observatory are 
also displayed (Abraham et al., 2008). 



on their observability. A non-exhaustive compilation of the 
UHE photon flux for some of the models presented above is 
shown on Fig. IT] The bottom panel shows the integrated fluxes 
compared to the current upper limits from the Pierre Auger 
Observatory (Abraham et al., 2008; see also Abraham et al., 
2009). The maximum energy was increased to 10^' eV for the 
mixed composition (the cosmic ray output then overshoots the 
data above 3 10'^ eV) and the late transition models to obtain 
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more optimistic fluxes than in Fig. |4||6] The three other mod- 
els (the second late-transition model with E,„„^ - 10^"^eV, 
the dip model and the pure iron model) are displayed with the 
default case E^ax = 10^'^^ eV. Evidently, the predictions dis- 
played here are far below the experimental upper limits. This 
means that, although these upper limits severely constrained 
the contribution of most (if not all) Top-Down models, UHE 
cosmogeiuc photons would be probably quite difficult to de- 
tect anyway (unless in the case of very optimistic models, e.g. 
a late-transition model with E„,ax > lO^^'^eV). These limits are 
expected to decrease by at least one order of magnitude dur- 
ing the Pierre Auger Observatory operating time (Risse and 
Homala, 2007). In this context the three most optimistic cases 
displayed in Fig. IT] corresponding either to very high values of 
the maximum energy or very hard spectral indices, might be 
constrained. 

It could be argued that a continuous source distribution is 
not a realistic modeling for the local universe. For a FR-II type 
sources, for instance, it is well known that this type of power- 
ful sources is extremely rare in the local universe. The cosmo- 
genic UHE photon flux would then be a lot lower than those 
displayed in Fig.|4]|6]and the observation of cosmic rays above 
3-4x10'^ eV would be difficult to explain (owing to the cos- 
mic ray horizon), while the low-energy photon flux and UHE 
neutrinos would basically remain unchanged. Furthermore, lo- 
cal over-densities (see for instance Gelmini et al., 2007; Taylor 
and Aharonian, 2010) in the source distribution as well as a 
nearby source with a high contribution to the total cosmic ray 
flux (see Taylor et al., 2009 and the next section) can strongly 
afffect the expected UHE photon flux. We study the impact of 
the contribution of single sources in section 6. 

5.3. Consequences of the UHE cosmic ray 
composition 

The discussion on the effect of composition on the expected 
flux of secondaries is complicated because the gap between the 
predicted fluxes of UHE neutrinos depends on the spectral in- 
dex required to fit the spectrum for a given maximum energy 
per unit charge, which itself depends on the source evolution 
model assumed. We have shown, however, that UHE photon or 
neutrino fluxes are quite similar for the pure iron and pure pro- 
ton cases. In any case, as soon as the maximum energy at the 
sources is above the pion production threshold (and the source 
luminosity evolution is non-negligible to provide a hard enough 
spectral index for a pure iron composition) the discussion of the 
detectability of UHE secondaries does not severely depend on 
composition. 

The Pierre Auger Observatory recently reported the largest- 
statistics composition studies above 10'^ eV (Abraham et al., 
2010) based on the energy evolution of the maximum of lon- 
gitudinal development of air shower (X,nax) and its spread. 
Obviously, the interpretation of this result in terms of composi- 
tion of the UHECRs must be considered with care because the 
hadronic physics taking place at first stages of the shower de- 
velopment is currently not well understood. However, the X^ax 
energy evolution behaves as if the composition were gradu- 
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ally becoming heavier with energy, from a light composition 
around the ankle to a much heavier composition above a few 
10''' eV. As argued in Aflard et al. (2008) and Allard (2009), 
such an evolution of the composition is difficult to justify above 
10'^ eV if all the different species present in the source compo- 
sition are accelerated to the highest energies (above 10^° eV). 
In this case, one would expect a composition that becomes 
lighter (for instance in the case of the galactic mixed com- 
position, see Allard et al., 2007) or more or less steady (low 
abundances of He and intermediate nuclei compared to heav- 
ier nuclei). The most likely explanation of the observed trend 
involves models where the maximum energy per unit charge 
is limited around Z x 10'^ eV or below, i.e. models assuming 
that the proton component is not accelerated at the highest en- 
ergies. A mixed composition model (enriched in iron to better 
match the experimental spectrum, Allard et al., 2008; Allard, 
2009) with E,nax = Z X 10''' eV (hereafter referred ti as the 
low-Einax model) was then proposed as a possible way to re- 
produce experimental spectra and the high-energy composition 
trend. Secondary fluxes for this specific model are presented in 
Fig. [8] Evidently, UHE photon and neutrino fluxes are depleted 
owing to the lack of cosmic rays accelerated above the pion 
production threshold and should be well below the detection 
threshold for all current experiments. 

In this version of the model, however, all sources are as- 
sumed to have the same maximum energy, which means none 
of them is able to accelerate cosmic rays above the pion pro- 
duction threshold. However, one could argue that a dispersion 
of the maximum energy at the sources (which was proposed by 
Kachelriess and Semikoz (2006) for instance, to make the dip 
model compatible with Fermi acceleration motivated spectral 
indices) is possible (if not likely). This point is usually illus- 
trated by the famous Hillas criterion (Hillas, 1984), stating that 
the Larmor radius of the accelerated particles cannot exceed 
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the size of the source (rL(E) < R,). One can simply relate the 
expected maximum energy at the sources with their magnetic 
luminosity. The simplest estimates (Achterberg, 2002) gives 
Emax ~ 2.5 X lO^^ZySsFs X (LB/10'*<'erg s-')'/^ ^ where yS^ and T, 
are the speed and the Lorentz factor of the shock (a much more 
detailed discussion can be found in Lemoine and Waxman, 
2009|j From this argument it is often stated that very pow- 
erful sources such as FR-Il galaxies or GRBs are prime candi- 
dates for the acceleration of cosmic ray protons above 10^° eV. 
This simple luminosity requirement is much looser for nuclei 
(especially with large Z), so one could easily conclude that a 
greater number of sources should be able to accelerate nuclei 
at the highest energies and that this component might have a 
dominant contribution to the total UHECR luminosity. On the 
contrary, powerful accelerators might be rare and outnumbered 
in the local universe or even absent from the UHECR horizon 
above a few 10'^ eV. Even if the current composition trend 
shown by the Pierre Auger Observatory is confirmed in the fu- 
ture, one cannot exclude, however, that powerful UHE proton 
accelerators might be able to contribute at the level of a few per- 
cent or a few tens of percent to the cosmic ray flux at 10'^ eV. 
Figures |4]|6]clearly show that, if these powerful accelerators are 
strongly evolving, like in the FR-II evolution model, a contribu- 
tion of 10% to the total flux at lO'** eV would be enough to pro- 
duce a diffuse neutrino flux detectable by IceCube and Auger, 
which would be a direct signature of the subdominant contribu- 
tion of these powerful accelerators^ This scenario would pro- 
duce a diffuse UHE neutrino flux as high as the expectations for 
the SFR source evolution scenario in the previous section. As a 
conclusion for this paragraph, even for models with low maxi- 
mum energies per unit charge (whatever the exact composition, 
because the absence of UHE secondaries is mainly caused by 
the low Emax and not by the source composition), UHE neutrino 
observations are still plausible and future more stringent upper 
limits will be critical to constrain the contribution of powerful 
accelerators. 



6. Photons and neutrinos as signatures of 
powerful cosmic ray accelerators 

Independently of the global composition at the highest ener- 
gies, the detection of individual sources accelerating protons 
above 10^" eV would represent a major step forward for un- 



^ Of course this is a necessary but insufficient condition because 
it is not deahng with energy losses during the acceleration process. 
However, adiabatic losses and synchrotron losses would provide larger 
maximum energy for nuclei than for protons (for synchrotron losses 
the scaling of En^x with A and Z depends on the energy evolution 
of the acceleration time). For losses caused by to interaction at the 
sources, this discussion strongly depends on the source environment 
(see for instance, Allard and Protheroe, 2009) 

■^ It is often said that the observation of UHE neutrino would inval- 
idate the X„,ax evolution reported by Pierre Auger as being a compo- 
sition feature and favor an interpretation based on a change of high- 
energy hadronic physics. Current and future constraints from LHC on 
high-energy hadronic phenomena taking place in air shower develop- 
ment (d'Enterria et al., 201 1) would be crucial to estimate how likely 
this alternative interpretation is. 
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Fig. 9. UHECR and UHE photon-integrated fluxes for a source 
located at 20 Mpc and diff'erent values of the maximum energy: 
2019.5^ lO^", lO^"-^ and 10^' eV. Top: a source spectral index 
/3 = 2.0 is considered. Bottom: j3 = 2.6, the integrated flux is 
normalized assuming that the source provides three cosmic ray 
events above 5 x 10'^ eV. 



derstanding the origin of UHECRs. Cosmic ray observatories 
are of course prime candidates for that purpose with the mea- 
surement of cosmic ray anival directions at the highest ener- 
gies. Current statistics is not large enough to firmly establish 
the detection of significant clustering caused by to the contri- 
bution of individual sources, but the next years should allow 
the Pierre Auger Observatory and JEM-EUSO (Medina-Tanco 
et al., 2009) to provide more detailed and extended statistics 
measurements. Above ~ 3 x 10^" eV heavy nuclei are expected 
to disappear from the cosmic ray composition (provided the 
closest source of UHECR is not much closer to the Earth than 
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Fig. 10. Same as Fig. |9| but we consider here different source 
distances (3.8, 20, 50, 70, 100 Mpc), the maximum energy is 
fixed at E^ax = lO^"'^ eV and the spectral index/? = 2.0. 



Centaurus A and the abundance of nuclei heavier than iron is 
low), owing to interactions with CMB photons (Allard et al., 
2007; 2008). Cosmic rays at these energies (if any) should be 
widely dominated by protons or very light fragments (i.e., parti- 
cles with very high rigidities). This energy range, which should 
be probed and constrained by JEM-EUSO, is then very promis- 
ing to observe cosmic rays that are almost guaranteed to point 
back close to their source. As a consequence of the cut-off of 
heavy nuclei, however, the cosmic ray flux might be very low 
if cosmic rays are dominated by heavy nuclei below 3 x 10^" 
eV. 

6.1. Cosmogenic UHE photons from single sources 

Cosmic ray observatories are also expected to spot individual 
sources by observing UHE photons. As mentioned above, UHE 
photon flux is strongly influenced by the distribution of sources 
in the local universe. Therefore, we studied, as in Taylor et al. 
(2009), the UHE photon flux of a source accelerating protons at 
UHE and located at various distances, with various hypotheses 
on the maximum energy. The example of a source located at 20 
Mpc and different values of the maximum energy is shown in 
Fig. l9j where we assumed that the source was providing three 
events per year above 5 x lO'^ eV (corresponding to a bit less 
than 10% of the flux seen by the Pierre Auger Observatory). 
Obviously, for a spectral index y6 = 2.0 and maximum energies 
above 10^°'^ eV the rate of expected events (these are incom- 
ing events and not events detected and identified as photons) 
is between ~0.4 and 1 per year above 10'^ eV and between 
0.08 and 0.2 above lO'^'^ eV, which is almost as high as the 
diffuse flux shown in Figs. HI and 17] for the dip model. The flux 



flux; the expected flux appears to drop significantly for D>50 
Mpc compared with closer sources owing to the cascading of 
the electromagnetic particles to lower energies. In contrast, for 
very nearby sources (for instance, we choose 3.8 Mpc, the dis- 
tance of Centaurus A) the flux above 10''^ eV is on the same 
order as in the 20 Mpc case, but harder because of the shorter 
propagation time that shifts interacting protons to higher ener- 
gies and prevent the cascades from developing too much toward 
lower energies. 

As commented above, if the source is close enough and 
has a significant contribution to the total UHECR flux (~ 10% 
above 5 x 10'^ eV) its contribution to the UHE photon flux 
can be as high as the whole diffuse photon fluxes calculated 
in Fig. 4]|6 above lO'^ eV, provided its maximum energy 
is around 10^"^ eV or above (a relatively hard spectral in- 
dex would also help). For a proton-dominated composition the 
presence of one or a few of this type of sources could greatly fa- 
cilitate the detection of UHE photons compared to the default 
case of a continuous source distribution that we studied. In a 
low Emax type scenario, the detection would be more difficult 
because a dominant contribution of sources with low values of 
the maximum energy would not provide a background of UHE 
photon that can be added to that of the single powerful acceler- 
ator and a higher contribution of the individual proton source(s) 
would be needed to compensate for the diffuse flux. A detection 
of UHE photons would be especially precious because it would 
prove an UHE accelerator above ~ 10^° eV per nucleon (which 
is currently not obvious) and it would be of special interest to 
compare their arrival directions with those of the highest en- 
ergy cosmic rays. 



6.2. Distant point sources of cosmogenic UHE 
neutrinos and GeV-TeV gamma-rays 

As mentioned above, most powerful accelerators could be lo- 
cated well beyond the cosmic ray horizon (this could be true 
even if the cosmic ray composition is proton-dominated at the 
highest energies). The only signature of the cosmic ray accel- 
eration taking place in these sources would be the observa- 
tion of high-energy neutrinos or sub-TeV gamma-rays that can 
propagate on cosmological distances, unlike UHE photons or 
cosmic rays. Studies on the subject have been most recently 
undertaken by Essey et al. (2010) and Ahlers and Salvado 
(2010) for nuclei. One potential problem of this type of de- 
tection is that distant sources have to be extremely powerful to 
provide detectable fluxes (see below) at least for cosmogenic 
neutrino. As an illustration, the example of a source at 1 Gpc 



with a cosmic ray luminosity of 10 erg s 



-v47 



' between 10'^ and 
for a pure proton 
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Emax = Z X lO^'^'^ eV is displayed in Fig 
and a pure iron source composition. The top panel shows the 
contribution of this source to the total cosmic ray spectrum as- 
suming a completely rectilinear propagation to ensure the very 
high luminosity assumed does not overshoot the UHECR flux 
(however, large powers of this magnitude are known for pow- 
erful astrophysical objects in other wavelengths). Obviously, 
the source contributes ~ 15%at~ 10'^ in the pure proton 
shows the influence of the source distance on the UHE photon case, before a cut-off caused by interaction with CMB pho- 



becomes obviously lower for softer spectral indices. Fig. 10 
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Fig. 11. Top: Cosmic ray spectrum at the Earth for a source 
located at 1 Gpc emitting a cosmic ray luminosity of 10"*^ 
ergs"' between 10'^ eV and Emax = Z X 10^°'^ eV assuming 
a pure proton (blue solid line) and a pure iron composition (red 
dashed line) compared with the Pierre Auger Observatory spec- 
trum. Bottom: Cosmogenic neutrinos (dashed lines) and pho- 
tons (solid lines) for the pure proton (blue) and the pure iron 
(red) cases corresponding to the same cosmic ray source. 



tons, whereas the iron source contributes at the level of 30%, 
their existence is then not constrained by the UHECR cosmic 
ray spectrum. 

The corresponding cosmogenic neutrinos and photons are 
displayed in the bottom panel of Fig. [TT] For such a high lu- 
minosity the neutrino flux can reach 10^ GeV cm"^ s"' at the 
peak around 3 10"* eV. A flux this high might be detected in 
the future depending on IceCube (Abbasi et al., 2009) or the 



expected Km3Net sensitivity to point sources at very high en- 
ergies. For lower energies the neutrino flux decreases rapidly 
and the detection of such a signal in the 1-100 PeV range would 
require a luminosity much higher (see Essey et al., 2010), es- 
pecially when using the Kneiske et al. (2004) estimate of the 
IR/Opt/UV backgrounds. The detection of a UHE neutrino sig- 
nal would provide an important non-ambiguous proof of the 
existence of very powerful UHE accelerators {i.e. acceleration 
above the pion production) even though the similar flux ob- 
tained for a pure iron composition, assuming the same lumi- 
nosity and same maximum energy per unit charge, shows that 
these fluxes are not sensitive to the source composition (at least 
not for hard spectral indices, see discussion above). Moreover, 
these neutrinos are produced in close vicinity to the source (at 
this redshift mostly within 10 Mpc) and are then very likely 
to form a point-like neutrino source for neutrino observato- 
ries. This is expected to remain true for beamed sources where 
the luminosity mentioned above would become the isotropic 
equivalent luminosity. 

Photons in the GeV-TeV range for single individual sources 
of UHECRs have akeady been considered in Ferrigno et al. 
(2005), Armengaud et al. (2006) and more recently Essey et 
al. (2010) and Ahlers and Salvado (2010). In this energy range 
the cosmogenic photon flux displayed in the bottom panel of 



Fig. 1 1 is very high and of course well within the range of ex- 
isting and future gamma-ray observatories. But ~ 75% of this 
flux are provided by the pair production mechanism, then most 
of the gamma-ray flux is produced far from the source (unless 
a very strong magnetic field is able to confine protons between 
10'^ eV and a few 10'^ eV). The sole flux caused by pion pro- 
duction remains ample however and far above the gamma-ray 
observatories sensitivity for point-source detection. As pointed 
out in Gabici and Aharonian (2007) (see also Gabici, 201 1) this 
flux is yet likely to be widely spread and isotropized as soon as 
the extragalactic field is not much lower than 10 '^ G. For mag- 
netic fields between 10"'^ and 10"^ G (above this value, syn- 
chrotron losses influence the cascade development, see below, 
and Gabici and Aharonian ,2005, 2007) as discussed in Gabici 
and Aharonian (2007) and assuming an isotropic cosmic ray 
emission, an electromagnetic halo of typical radius ~ 20 Mpc, 
including most of the flux produced by the pion production 
mechanism, would form. It would spread the source image over 
an angular size on the order of ~1 degree or more (depending 
on the precise value of the field, which is needed to estimate the 
3D development of the electromagnetic cascade) for a source 
located at 1 Gpc. Then, the angular size of the source would 
exceed the typical point spread function of a few arc-minutes, 
which should be chosen for the future CTA gamma-ray exper- 
iment, but the high expected integrated flux might still be de- 
tectable. 

For greater values of the magnetic field, especially in the 
source environment, Gabici and Aharonian (2005) pointed out 
the possible detection of the synchrotron photons emitted dur- 
ing the first stages of the cascade development. This possibil- 
ity was also recently studied in Kotera et al. (2011). In both 
studies, a source at 1 Gpc of luminosity lO'^^ergs"' above lO'^ 
eV is found to possibly yield detectable fluxes for Fermi or 
CTA (for instance, a flux of 10 '^GeVcm^^ s ' spread over 
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ity of lO'*'' X (D/1000Mpc)2ergs"'. The top panel shows the 
photons and neutrino fluxes, the bottom panel the spectra com- 
pared with the Pierre Auger Observatory data. Only the pure 
proton case is considered. 



~ 0.2 deg). These synchrotron images are promising for the 
detection of distant UHECR accelerators and have the advan- 
tage that they are a signature of the pion production mechanism 
(see discussions in Gabici and Aharonian, 2005 and Kotera et 
al. 2011), unlike cascades in the case of distant sources (see 
below for the local universe). However, their signal is weaker 
in the TeV range than in the GeV, although their detection with 
CTA should still be possible. However, if the matter density is 
quite high in the magnetized environment (see the example of 



a galaxy cluster in Kotera et al., 2009), the synchrotron signal 
could be screened and overwhelmed by hadronic interactions 
of the confined low-energy cosmic rays. 

Interestingly, GeV- TeV cascades were recently proposed by 
Essey et al. (2010) (and were recently re-examined by Alhers 
and Salvado, 2011) as a possible interpretation of TeV sig- 
nal observed by HESS of the AGN 1ES0229+200 (z=0.14). 
Assuming a point-source cascade signal, Essey et al. estimate 
the implied source luminosity in cosmic ray protons above lO'^ 
eV to be between ~ 10"*^ and 10^'' ergs"' (depending on the 
maximum energy assumed). If confirmed, this interpretation 
would not be a signature of the acceleration of particles above 
the pion production threshold (because this signal could be pro- 
duced by the sole pair production mechanism), but would in- 
volve either extremely weak magnetic fields (< 10 '"^ G) on 
very large scales or an even greater source luminosities (to al- 
low the same flux to be within HESS point spread function). 

One could argue that because the neutrino flux forms 
quickly, similar fluxes could be found by assuming weaker 
source luminosity Ld - 10'*'' x (D/1000)^ergs"'. Neutrino 
and photon fluxes for sources located at distances D=5, 10, 20, 
50, 100, 200, 500, 1000 Mpc assuming a luminosity Lq above 
10'^ eV, are shown in Fig. 
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For neutrinos only, fluxes for 10 
and 1000 Mpc are shown for clarity. Neutrino fluxes become 
completely similar to those at 1 Gpc above ~ 20 Mpc, slightly 
shifted to higher energies because of the lower redshifting of 
the neutrino and the higher energy threshold implied by the 
smaller distance of the source. These fluxes requiring lower lu- 
minosity could then be detected if neutrino observatories were 
able to reach a sensitivity around 10"'^ ergs"' above lO'^ eV. 
However, as can be seen in the bottom panel of Fig. [12] unlike 
for a source at 1 Gpc, for which energy losses keep the prop- 
agated spectrum safely below the experimental data, a source 
located at D< 500 Mpc would overshoot the cosmic ray spec- 
trum for luminosities on the order of L/Jj For a reasonable 
contribution to the cosmic ray spectrum, say on the order of 
10% around 5 x 10''^ eV, one would have to assume a lumi- 
nosity ~ 20 times lower, e.g. around ~ 5 x 10^' ergs"' above 
10'^ eV (assuming/? = 2.0). Then, the flux would drop around 
5 X 10"'° GeV cm"^ s"', which should be a lot more difficult to 
spot. Yet, reaching this sensitivity would make neutrino obser- 
vatories extremely constraining for the potential accelerators in 
the local universe. For the time being, the detection of an ex- 
tremely powerful and distant source (which might be the prime 
candidate for proton acceleration above 10^" eV, see above) 
seems to be a more realistic goal. 

The corresponding cosmogenic photon fluxes (UHE and 

The 
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GeV- TeV) are also displayed in the top panel of Fig. 
GeV- TeV flux becomes lower as the source distance decreases. 
This is because the contribution of the pair production mech- 
anism decreases. The higher contribution of the pion decay 

■* For the farthest distances it could be argued that the limited life- 
time of the source (typically a few million or tens of million years for 
AGNs) and a hypothetical extragalactic magnetic field might cause a 
dilution of the flux during the propagation. This is quite unlikely for 
sources located at less than 100 Mpc, however, and our discussion 
cannot include bursting sources. 
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increases the average energy of secondaries that initiate the 
cascade and consequently the electromagnetic cascades take 
longer to develop down to the TeV range (magnetic field ef- 
fects would of course make this discussion more complicated) 
which in turn makes the spectra harder, but the fluxes lower. 
As mentioned before, to obtain the sources that contribute to a 
level of ~ 10% of the flux around 5 x 10'^ eV, one needs fluxes 
lower by a factor of ~ 20. In the last case, they are likely to re- 
main above the future CTA point source sensitivity. However, 
as mentioned above, (and following the discussion of Gabici 
and Aharonian, 2007), the effect of a potential extragalactic 
magnetic field on the angular size of the gamma-ray image of 



the source is expected to be much worse for nearby sources, 
making the detection of their signal quite hypothetical. 

Finally, ignoring the extragalactic magnetic fields, we con- 
clude with a discussion on the general impact of the compo- 



sition and the maximum energy. Fig. 13 shows the cosmic ray 
output (top panel) and the cosmogenic photon flux assuming a 
source located at 20 Mpc and a luminosity L = 2 x 10"*^ erg s"' 
between 10'^ eV and E^j^x. We studied two source composi- 
tions, namely pure proton and pure iron, and two maximum en- 
ergies, Zx 10^°'^ eV and Zx 10'^ eV. As mentioned above, with 
this luminosity the contribution to the total cosmic ray spec- 
trum becomes plausible (between ~ 10 and 20% at 5 x 10'^ eV 
depending on the model). For the highest maximum energy the 
pure proton and pure iron fluxes are once again very close at 
low energies. For the low E^ax sources (Emax = Z x 10''^ eV, 
well below the pion production threshold, the proton flux drops 
by a factor of ~ 5 compared to the high E^ax case. The drop 
in the flux is, in this case, limited because the pair production 
mechanism, although less efficient for nearby sources, emits 
lower energy secondaries that cascade more rapidly to TeV en- 
ergies, whereas the high-energy secondaries resulting from the 
pion production mechanism are more likely to remain at high 
energies if they are emitted close to the observer. In the low 
Emax scenario for the pure iron composition, the fluxes are even 
higher and closer to high Emax scenarios. Indeed, besides sec- 
ondary protons, which are efficiently emitted by iron nuclei 
close to the maximum energy, the flux is produced quite effi- 
ciently by the pair production mechanism above ~ 6 10'^ eV 
with a loss length ~ Z^/A times shorter than for protons at the 
same Lorentz factor (see Allard et al., 2006, for the contribu- 
tion of the different energy loss processes on iron nuclei and 
discussion in Alhers and Salvado, 2011). The precise differ- 
ence between the different compositions or maximum energies 
of course depends on the assumed spectral index, but these ex- 
amples show that besides the probable angular spread that is 
likely to prevent the detection of these TeV images, GeV-TeV 
fluxes suffer to some extent from a degeneracy on physical pa- 
rameters such as the composition (which should also apply to 
neutrino sources) or maximum energy, the latter preventing in 
principle the use of these fluxes as a signature of sources ac- 
celerating cosmic rays above the pion production threshold, 
~ 10^° eV per nucleon (unless a neutrino, UHE photon, of a 
UHE cosmic ray counterpart is also found). On the other hand, 
a positive detection would bring unprecedented constraints on 
the extragalactic magnetic fields in the local universe and cer- 
tainly allow a clear identification of the source. 

7. Conclusion 

We have considered the production of cosmogenic secondary 
photons and neutrinos during the extragalactic propagation of 
UHE protons and nuclei. We discussed the constraints obtained 
from the Fermi observations of the diffuse gamma-ray flux. We 
found that, because the possibility of a high contribution of 
UHECR to the diffuse gamma-ray flux is cutTently not ruled 
out, significant UHECR neutrino fluxes observable by present 
observatories such as IceCube or the Pierre Auger Observatory 
can still be expected and would be of prime interest for the 
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field. The UHE cosmogenic photons could also be expected for 
optimistic astrophysical assumptions. These fluxes are not con- 
strained by the experimental diffiise gamma-ray flux because 
they show little dependance on the evolution of the source lumi- 
nosity evolution unlike UHE neutrinos. We also discussed the 
influence of the composition evolution at the highest energies. 
We pointed out that scenarios for which the maximum energy 
at the source is limited (well below lO^'' eV per nucleon) are ex- 
pected to have their UHECR neutrino flux strongly suppressed 
above 10'^ eV. However, we found that the likely variation of 
the maximum from source to source could yield an observable 
neutrino flux, if rare powerful accelerators (e.g, source able to 
accelerate particles above 10^" eV per nucleon) show a strong 
cosmological evolution of their luminosity and contribute at the 
level of ~ 10% to the cosmic ray flux around 10''^ eV. 



We also considered the possibility of observing cosmo- 
genic secondaries from individual sources to constrain the ori- 
gin of highest energy particles. Distant and very powerful 
sources were found to be prime candidates, and neutrino flux 
along with synchrotron images (Gabici and Aharonian, 2005) 
would be unambiguous signatures of the acceleration of par- 
ticles above the threshold of pion production. These obser- 
vations would not strongly constrain the source composition, 
though. But observations of sources in the local universe would 
require, with the current experimental capabilities, too high 
contributions of these individual sources to the UHE cosmic 
ray spectrum. This can be alleviated, however, with improved 
sensitivities (it would require point source sensitivities below 
10"^GeVcm"^ s"' above 10'^ eV for neutrino observatories, 
for instance), which would at the same time decrease the lu- 
minosity requirements for distant sources. GeV-TeV photons 
from electromagnetic cascade, taking advantage of the upcom- 
ing very sensitive CTA observatory, could be detected even in 
the local universe by assuming a straight line propagation of the 
electromagnetic cascades. However, their detectability strongly 
depends on the strength of the extragalactic magnetic field and 
sources in the local universe are likely to become very hard 
to detect if the latter exceeds ~ 10~'^ G in the local universe. 
For a positive detection strong constraints on the extragalactic 
magnetic field could of course be obtained as well as a precise 
location of the source - a critical milestone in high-energy as- 
trophysics. However, these cascades are not a signature of the 
acceleration of cosmic rays above the threshold of pion produc- 
tion or of the source composition. Consequently, counterparts 
(primarily from cosmic ray observatories) would be needed to 
fully understand the signal. 



We conclude that the observation of cosmogenic sec- 
ondaries either by their diffuse flux or individual sources 
can be expected from current or next-generation instruments. 
Importantly, they would represent invaluable additions to ob- 
servations of the UHECR sky made by the Pierre Auger 
Observatory, the Telescope Array and the future JEM-EUSO 
in their quest for solving the long-standing UHECR problem 
and identifying its astrophysical source. 
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